Perturbation Methods and First Order Partial 

cn" : Differential Equations 

o 

David Holcman * Ivan Kupka^ 

o 

Q - 

Q ; November 26/2002 



-i— > 

B 

> 



Abstract 



> 

m 

I In this paper, we give explicit estimates that insure the existence of solutions 

■ for first order partial differential operators on compact manifolds, using a viscosity 

method. In the linear case, an explicit integral formula can be found, using the 
characteristics curves. The solution is given explicitly on the critical points and the 
limit cycles of the vector field of the first order term of the operator. In the nonlin- 
O ear case, a generalization of the Weitzenbock formula provides pointwise estimates 



that insure the existence of a solution, but the uniqueness question is left open. 
Nevertheless we prove that uniqueness is stable under a C 1 perturbation. Finally, 
we give some examples where the solution fails to exist globally, justifying the need 
to impose conditions that warrant global existence. The last result reveals that the 
zero order term in the first order operator is necessary to obtain generically bounded 



^ ■ solutions. 
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1 Introduction 



In this paper, we consider a compact orientable Riemannian manifold (V,g), a smooth 
vector field b on V or a parametrized vector field b:IRxV — >TV, the tangent bundle of 
V, and a smooth positive function c:V — >IR or a parametrized positive smooth function 
c:MxV — >IKL A g is the negative Laplacian: 

A g = —div g grad g 

where div 9 denotes the divergence operator with respect to the volume form associated 
to the metric g and grad ff the gradient with respect to g. 

We shall study the limit of the solutions u e as e tends to through positive values of 
the linear equations: 

eA g u e + < b,gradu £ > g +cu t = f on V (1) 
where / is a given continuous function on the manifold V or of the non linear equations: 

eA g u t (x)+ < b(u £ (x) , x) , gradu e (x) > g +c(u £ (x), x)u £ (x) = 0, x e V (2) 

Heuristically, the limits if they exists, "should" be solutions of the first order partial 
differential equation: 

< b, gradu > g +cu = f on V (3) 
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in the linear case and: 



< b(u(x) , x) , gradu(x) > g +c(u(x), x)u(x) = 0,x e V (4) 

in the nonlinear case. The linear case is obviously a particular case of the non linear one. 

Historically, Cauchy devised a powerful method to find local solutions of the Cauchy 
or initial value problem for first order partial differential equations (linear or non linear) 
using the characteristic curves which were the solutions of an ordinary differential system 
in the first jet bundle, called the characteristic system (see [3]). The initial data were 
given on a hypersurface in the bundle of the first jets of functions on the space of the 
partial differential equation and then propagated along the characteristic curves from 
the hypersurface. In the case of equation (4) , using the natural coordinate system 
(x 1 ,...,x m ,u,pi,...,p m ), on the jet space, the characteristic system for equation(4) can be 
expressed as follows: 

f = (5) 
k=i 

dp n nt^ ( db h , db k , A , dc dc , . 

-m = - 1^ [ dx^^ x) + Pn ^ u ' x n^- p»( u fa( u > x ^ + c ^ x ^ - x ^ 

k=l V ' 

where b k (u,x),l< k < m, are the components of the field b. 

The problem of solving (4) is not an initial value problem but a Dirichlet problem. 
It is not well behaved because the differential operator appearing in the left hand-side of 
equation(4) is not elliptic. Nevertheless the Cauchy characteristics are still very useful for 
the study of equation(4). Indeed in the linear case one can give an explicit formula for 
the solution u of equation(3) : 

/o 
f(x P (s))e-^ c{xp{r))dT (6) 
-oo 

see formula(lO) in theorem(2) below. Actually this formula can be used to prove the exis- 
tence and properties of classical (i.e. at least Lipschitz continuous) solutions to equation(3) 
due to the fact that the equation(5) can be solved without the knowledge of the solution 
u. Actually, provided that minc>0, the formula (6) defines a function u on V but this 

function is, in general, not even continuous. In the non linear case, the formula(6) is 
replaced by the integral equation in u: 

,o 

u(x) — / K(t,x)dt, x G V, where 

J — oo 

K(t,x) = f(v?"(t, x)) exp[- c(m(v? u (s, x)), y? u (s, x))ds] 
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see equation(33) below. This integral equation can be solved using the Picard fixed 
point theorem and give us a Lipschitz continuous solutions under assumptions which are 
stronger than the ones needed to prove the existence using the viscosity method developed 
below. 

The standard viscosity method is the approach using an elliptic partial differential 
equation to find solutions for a first order partial differential equation. Other methods 
such as the one mentioned in the last paragraph, exist but the advantage of the viscosity 
method is that, explicit geometrical conditions and explicit constants in inequalities are 
derived. Under these explicit conditions, the existence and local uniqueness of solutions 
are established. We prove that in certain types of first order linear partial differential 
equation, the solution is unique and a representation of this solution can be found using 
the trajectories of a vector field. 

First order partial differential equations similar to equation 4 and 3 were considered in 
relation to the KAM theory by Kuksin in [16]. Here V is a m-dimensional torus endowed 
with the flat metric b is a translation-invariant purely imaginary vector field, c must 
satisfies the condition (dx is the Haar measure) 



c(x)dx 



v 



» max{osc(c),l}. (7) 



Under some incommensurability assumptions on the coefficients of the field in the standard 
basis of the lie algebra of the torus and the condition 7, Kuksin proves that equation 3 
has a unique analytic solution for any analytic f. Moreover sup-norm estimates of the 
solution are given in terms of the sup-norm of f. The method of proof is based on the 
theory of similar equations on the circle S 1 . Because we are on a torus Fourier series can 
be used to represent the functions and this is how the estimates are derived. 

Forni in [6] studies equation 3 with c=0 on compact surfaces endowed with a Rieman- 
nian metric. He assumes that the vector field b preserves the volume measure associated 
to the metric and that its singularities are nondegenerate saddle points only. When f is 
sufficiently regular and its mean value vanishes, there exists a unique solution( up to a 
constant of course). 

In the following paper, there are no restrictions on the underlying manifold V (except 
compactness), b is a general Morse-Smale vector field. Hence in particular it is not volume 
preserving. Also in our case c cannot be taken to be 0. It has to satisfy condition min 

c> b (defined in the notations below) which is more stringent than 7. This condition 
means that c must be larger than the variation of b. In the case considered by Kuksin, 
this variation is because b is self-parallel. 

Some of the results exposed here were announced in [9]. Let us finally mention that 
when the function f in the right hand side of equation (2) is proportional to u, we get an 
eigenvalue problem. The analysis of this problem has been done in the general case in 
[10, 11]. 

The paper is organized as follow. In the first section of the paper we study the limits 
as e goes to of the solutions of equation (1) to prove the existence and uniqueness of so- 
lutions for the linear equation (3). Then we prove the uniqueness of the solutions. In the 
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second section we take up the existence and uniqueness of solutions of the nonlinear el- 
liptic equation (2). The estimates needed to achieve this are provided by the Weitzenbock 
formula and its generalization to covariant tensors of valence two. Then as in the first 
section we examine the limits of solutions of equation (2) as e goes to and prove the ex- 
istence of solutions for the nonlinear first order equation (4). In the third and last section, 
we give various examples in the nonlinear case where a solution can exists in an open set 
but not globally. When the field is ergodic, we give a simple result about the behavior 
of the solution on the manifold. Finally in order to understand how relevant is the zero 
order term, we consider the case where it is zero. In that case, the viscosity method does 
not converge to any bounded solutions. Generically, a first order linear partial differential 
equation with no zero order term cannot have bounded solutions. 



1.1 Notations 



Throughout this paper we shall use the following notations: 
l)Metric and tensors: 



d g = distance function defined on V by the metric 
|| * ||, < *, * >g— norm and scalar product associated to g 
volg = volume measure associated to g 
grad = gradient operator associated to g 
divg = divergence operator associated to g 
V = Levi-Civita covariant derivative associated to g 
R = curvature tensor of 
9(b) := Lie derivative operator associated to the vector field b 
9{b)u = du(b) 

For any coordinate system x 1 ,.---,*" 1 : O — >R on V : 

g = gijdx 1 <g> dx 3 
g %3 = inverse of the matrix g i ■ 



vol g = ^Jdet(g ij )dx 1 ...dx r 



gradu = g J 



dx'j 

diVgX = V i X i 



d 

R = R-l dx 1 <g> dx 3 (g)dx k (g)—r 
•j* q x i 
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If T is a p-covariant,q-contravariant vector field: 

T = t*;;;;*" •-%/./•' <g> ... <g> g^ +1 <g> e h <g> ... <g> e i9 

VT = V^'^-^W 1 <g) ... <g) cfe ip+1 <g> <g> ... <g> e jq 
where (ei,...,e m ) is the frame field associated to the coordinates x 1 ,.--,* 7 ™, on O. 

2)Norms: g induces a scalar product function and a norm function on any tensor 
bundle on V. Let t, r be two tensors of the same type, at the same point of V. 

< t,r > g = scalar product of t and r 
| \t\\ g = norm of t 

To any tensor field T on V is associated the function: 
xG V— >\\r \\ g E [0,+oo[. 

Halloo = max||T (x)|| 
'|^~Hl°° = ess-sup||T (x)|] g 



r \\tf ={ W r (x)\\ q vol g (dx 



V 



W\W = Halloo + I|vt i 

|7"||i4/i,oo = llT'll^oo + ||V7" 



oo 

I Il c 



\nm = m\ni 2 + \\vr\\i 2 



3)Constants: 



6 = sup{< V x b,X > g \X e TV, \ \X\\ g = 1} 

, . db . . . . db . . 

P = hrr °° = max hrr 9 
OX Kxv d\ 

en = inf c(A, x) > 

AeiR,xeV 



To = |]a/||oo + - 



c 



_ 1 1 9C . |./ ;, x 



<9A c 

r = max{ | Ricc(u, uj) | ,| w G T*V,||u;|| s = 1} 



4)Lipschitz properties 
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A function f:V — >R, will be called Lipschitz continuous if 



\f(x)-f(y)\ ^ 
sup — — < +00 

x>yeV d g {x,y) 

The sup is called the Lipschitz constant of f and will be denoted by Lip(f). Any upper 
bound of the Lipschitz constant is called a Lipschitz bound of f. Finally we denote 



C 0,1 (V) = Banach space of all lipschitz continuous functions on V with the norm 
||/||co,i = max|/| +Lip(f) 

2 The linear case: main theorem 



In this section, c:V — >R, will denote a smooth positive function and f:V — >R, a Lipschitz 
continuous one. Let us recall that by the theory of second order elliptic equations ([2]), 
the following equation: 

eA g u e + < b, gradu £ > g +cu e = f on V (8) 

has a unique solution u e which is of any class C 2 ' Q , for a G [0, 1[. If f is of class C k , 
ke NUoo, u e is of class C k+2 . We shall study the limit when e converges to zero of the 
solution u e oi the equation(8). 

Theorem 1 1. On a compact manifold, consider a smooth vector field b, a smooth 
positive function c satisfying mine = Cq > b and a Lipschitz continuous function 

f. Under these assumptions, there exists a unique Lipschitz- continuous function 
u: V — >~R, solution of the first order partial differential equation : 

< b, gradu > g +cu = f on V (9) 

in the weak sense. 

2. For any P& V, if xpM — > V, denotes the trajectory of b such that xp(0)=P,we have 
the formula: 

/o 
f(x P (s))e'^< Xp ^ ds dt (10) 
■00 

3. The solution u £ of the equation eA g u t + < b, Vu e > +cu e = f tends to u in the sup 
norm topology as e tends to zero. 

4- If f is of class C 1 this solution is also C 1 . 

5. At a singular point P of b, u(P)=^^. 
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6. If p:M — > V n is a periodic trajectory of b with minimal period T, for all tE R 

u(p(t)) = f /««»gf* + T -i^ jf M a« t, 

(11) 

where C(t)—exp[—J Q c(p(s))ds\. 

Proof : First we will prove that any sub-sequences of u t contains a converging sub- 
sequence whose limit u in the weak L 2 topology, is a Lipschitz continuous function sat- 
isfying equation (9). In order to show this we need to establish some a priory estimates 
which follow easily from the maximum principle. The following is easy: 

min / max / 

— ^- <u e <— v — (12) 
max c mm c 

v v 

There exists a sub-sequence of u e which converges weakly to a function u and for all 

e / u e A g (j) — / div((f)b)u e + / <pcu e = j 4>f 
Jv Jv Jv Jv 

u satisfies for all G Hi, 

div{4>b)u + / 4>cu — I (j)f 
Jv Jv 



In these equations the integrations are with respect to the measure defined by the metric 
g on V. 

Now we will prove that | \du e \ ^ is also bounded. Let us note here that using a stronger 
lower bound on c, one could prove that 1 1 S7du e \ \ g is also bounded in which would imply 
differentiability of the limit. To estimate 1 1 rfit e 11^,, we shall consider the maximum of the 
function ||cfoi e || 2 , following a standard method (see [17, 13, 19] ). 

Starting with the Weitzenbock formula: 

-A ff ||dw e ||g + ||Vdu e ||g + Ric(du e , du e ) =< dA g u e , du e > g 
2 



Assuming that f is C 1 , an easy computation gives 



e ]-A g (\\du e \\ 2 g ) + e||Vc£u e ||g + eRicci(du e , du e ) + c||rf-u e || 2 + < V gra du f b, gradu e > £ 

1 



= — -^9(b)\\du e \\ 2 + < df,du e > g —u e < dc,du € > g . 
2 tf 

Because the Ricci curvature is bounded from below by r and u e is bounded from above 
, we obtain that at a maximum point P of the function ||dw e || 2 



(c - er - 6o)||d« e (P)||^ <|< df(P),du e (P) > g \ + \u £ (P) < dc(P),du e (P) > £ 
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Note that at a maximum point P of the function ||du £ ||^, A g (\\du t \\ 2 g )(P) > and 
(9(b)\\du\\ 2 g )(P) = 0. Hence we get using the relations (12): 



\W\\ 



wl = n^)ii fl < : fr - h (13) 

y c — er — o 



Now we drop the auxiliary assumption that fe C 1 and prove that for any Lipschitz con- 
tinuous f: 



\duA\ < " " 



c - er - o 

We use the following lemma (probably well known) proved in the appendix: 



Lemma 1 Let k:V — >M ; be a Lipschitz continuous function with M as Lipschitz bound. 
Then for any neighborhood it of k in CP(V), any e > 0, there exists a C°° function h:V— 
>R, contained in it and admitting M+e,as Lipschitz bound. 



Hence with the lemma, we can find a sequence {f n \n G N} of smooth functions on 
V, converging to f in C°(V), f n allowing ||d/||Loo + - as Lipschitz bound for each n. Let 
u e>n be the solution of equation (8) with second member f n . It follows from classical 
elliptic estimates that the sequence {u ei „|n G N} converges to u e in the C 2 topology. The 
inequality (13) implies that for each nG N: 



\du fn \\_ < 



1 1 dfn I I oo 4" 



CO 



Co - er - b 

for every e G [0, C() ~ b ° ]. From this it follows that u £ is Lipschitz continuous and that: 

| M/ || L ^ + MkJMk 



|diiJ|„ < 



c - er - b 



Because \\u e \\ g ,\\du e \\ g are uniformly bounded in e for e G [0, c °~ b ° ], any sequence contains 
a sequence converging uniformly and even in any Holder space C°' a , for a < 1, to a solution 
u of the equation 

< b(x) , gradu(x) > g +c(x)u(x) = f(x) almost everywhere on V 

Also ||du||L°o is finite and: 

\\df\\ L ~ + m il dcll °° 

\du\ r oo < : • 

c - o 
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2.1 proof of existence and uniqueness 



Now we prove that there is only one function u:V — > R, which is a weak solution of 
equation(9) and is Lipschitz continuous. We have just shown that such u's exist. Taking 
any one of them we will give an integral representation of the function u along the tra- 
jectories of b which will be used to compute the value of the function u at the stationary 
points and along the periodic trajectories of b. Since u is Lipschitz continuous, it is abso- 
lutely continuous along any C l curve and almost everywhere differentiable. If x:R — > V 
is any trajectory of b, the function: te R — > u(x(t)) is absolutely continuous and satisfies 
the equation: 

du(x(t)) 

: h c(x(t))u(x(t)) = f((x(t)) almost everywhere (14) 

at 

Hence the functiomtG R — > u(x(t)) is a genuine C 1 solution of equation(14) and we have 
for all t, ti, 

u(x(t)) =u{x{t 1 ))e-^ cix(s))ds + f f{x{s))e-S*« T » dt ds 
Because u is bounded and mine ^ Co > 0, we get, if we let ti go to -oo : 



l (x(t))= [ f(x( 8 ))e-X**V 

J — oo 



For any point P in V, if x P :R — >V denotes the trajectory of b passing through P at time 
0: 

/o 
f(xp(s))e-I' c(lp(T))(iT 
-oo 

This formula proves the uniqueness of u. It also proves the assertion 3) 



2.2 proof of regularity 

We need an elementary lemma: 

Lemma 2 1. If p is the modulus of an eigenvalue of the linearization of the flow at 
time one of b at a singular point of b, then \ logp \ < b . 

2. If p is the modulus of an eigenvalue of the monodromy of a periodic orbit of b of 
minimal period T, then \ logp\ < b T. 

3. If {4>t I t E R} denotes the flow of b and T(p t its tangent mapping, then \ T(p t \< e b °^ 
for all real t. 
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Proof of the lemma: The first statement is easy and left to the reader. Let us prove 3 
first. Denote by v any tangent vector to V. Let v : R — > TV, be the vector field : 



v(t) = T<Mv ) 

We have: 

V t v = V v b 
d\\v(t)\\l 



dt 

d\Ht)\\: 



= 2 < V v(t) b,v(t) > s 
<26d|K0IIJ 



dt 

From this it follows that if t> 0: 

\Ht)\\ 9 <\H0)\\ 9 e bot 
Reversing the time, we get for t < 0: 

\\v(t)\\ 9 <\H0)\\ 9 e boltl 

This last inequality implies 3). 

Let us prove 2). Let p:M — > V be a periodic trajectory of b of (minimal) period T and let 
A be an eigenvalue of the monodromy of p. There exists a tangent vector v GT p ( )V®C, 
v 7^ 0, such that T0 T (v o )=Av o and T0_ T (v o )=iv o . By the above: 

|A| <e boT 
1 



< e b ° T 
|A| " 



This two inequalities imply assertion 2). 



Now we resume the proof of assertion 4). Let {<f>t\t G M} denote the flow of b. 
(p(t,x) = (p t (x). 
For any xGV: 

u{x) = / K(t,x)dt (15) 



where K:R x V — >■ R, is the function: 

K(t,x) = /(0(t,a;))e-^ OcWs ' a;))ds 
Let K t (a;) =K(t,x).K is a C°° function. Its differential in x is: 

dK t (x) = e-ti^^WWt^))^ - J° dc{<j>{s,x))T x <i> 9 ds] 
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T x (f) t is the tangent mapping of the diffeomorphism </?<:V — >V, at x. The second integral 
is that of the curve: sG R — > dc(4>(s,x))T x (j) s GT* V, the cotangent space of V at x. Using 
the statement 3) of the Lemma(2), for t< 0: 

\\dK t (x)\\ < e^-^iWflU + 1*111011!] 

This estimate, uniform for all x GX and all t< 0, shows that u is continuously differentiable 
on V because of the assumption Co >t>o and that: 

,o 

du(x) — / dK(t,x)dt 



2.3 explicit formulas on the recurrent set 



Now we can specialize the formula(15) to different kinds of trajectories. Let P be a 
singular point of b. Then the curve x:R — > V, x(t) = P for all t is a trajectory and 

u(P) = u(x(0)) = f f{x{s))e- f f{P)e-S!«F>*r=W 

J -co J -co C{r ) 

Now assume that p:R — > V, is a periodic trajectory of (minimal) period T. Applying the 
general formula to p, we get 

u{p{t))= f f{p{s))e-5l<^ 

J —CO 

We choose te[0,T]. Then: 

/(p(*))^_v^ r nT f(p(s)) ds 



J -co C(S) 1 ^7-(n + l)T C(S) 

f- nT f(p(s)) dc _ c -(n +i) Cm f /(p(*)) 
J-(n+l)T C(s) J C(s) 

J -co C(s) 1 - e-cm J Q 



ds 
ds 
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Remark 1 : Let us define the following real number: 

B = mib (g)= inf max{| < Vxb,X > g \X G TV, \ X \— 1}, 

g g^Riem 

where Riem denotes the Riemannian structure on V. Lemma 2 shows that B is strictly 
positive (bigger than the minimum of the logp and ). For each function c such that 
infc>B , the formula(15) defines a function u c which is a regular solution of equation(9): 

indeed one chooses a Riemannian metric g on V such that infc > b(g) and applies what 

has just been proved. When infc= -Bo, the formula (15) still defines a function u c , but 

it is not even sure that it is continuous, even less that it is a solution of equation (9) in 
some sense. Hence the following question: what happens in this situation? What type of 
solutions, if any does the equation(9) has ? 

Remark 2 : Can one find an explicit expression for the limit of the solution u, on a 
general recurrent set? what about the case, when the recurent sets have an invariant 
measure? this will extend the formula obtained on the limit cycles. 

3 The nonlinear case 

In this section, we consider a parametrized vector field :A G R — >b A GVF(V), space of 
all vector fields on V. b:MxV — >TV will denote the mapping: (A, x) — >b(A, x) = b\(x).c: 
IRxV — >R, will denote a smooth function and c A :V — >K will be the function c x (x) = 
c(A, x). 

We shall study the following equation on V: 

< b u , gradu > +c u = f on V (16) 

where u:V — >M, is the unknown function, f:V — >R, is a given smooth function on V, b u 
is the vector field on V,b n (x)=b(u(x),x), c u the function c u (x)=c(u(x),x). 

We will impose some conditions on /, c and b in order to insure the existence of regular 
solutions. Let us note here that there are examples with discontinuous solutions. This 
has been discussed in [12]. To prove the existence of solutions for equation (16), we use 
a standard elliptic regularization techniques and proceed by successive approximations. 

The following two assumptions will be in force in the following sections: 

1. c — b — 7 > 0. This implies that c >0. 

2. (c -6 - 7 )2-4/ /3>0 

These assumptions say that the minimum of the function c must be large enough 
compared to the maximum dilation of the field b and are called conditions of hyperbolicity. 
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3.1 Existence of the solution for the elliptic PDE 



To prove the first part of the theorem, about the existence of a solution to the first order 
partial differential equation, we will built a sequence of solutions of some nonlinear elliptic 
partial differential equations with a small viscosity coefficient and later on, we will let this 
small parameter converges to zero. 

Now the assumptions (1-2) imply that there exists two strictly positive roots of the 
second order equation: f3X 2 — X(c — b — 7) + / = 0. The same property will be true 
for the equation: (3X 2 — X(c — b — er — 7) + fo = provided that 0< e < e, where 

E = r, (17) 

For e e [0, e[, we denote by R(e) > its smaller root. It is easy to see that R(e) is an 
increasing function of e. 



Estimates 

Proposition 1 On a compact Riemannian manifold (V,g), under the assumptions 1-2 
above, for e > 0, the sequence^ Uk\ k<E Z + } defined by the first term u = and the 
recurrence relation: eAuk+i+ < b Uk , gradiik+i > g +c Uk Uk+i = f, converges in the topology 
to the solution u £ of the equation: 

eA g u e + < b Ue ,gradu t > g +c Us u t = f 

Moreover for any e e]0, e[, 

IWL<B(e)- 

For any e in ]0,e[, ||c?u £ ||oo < R(so) for all e in ]0,£ [- 

Proof: The sequence Uk is well defined, using the results of the previous paragraph 
on the linear case. By the maximum principle: 

min / max / 

v ^ ^ v 
< u f < — ■ 



max c mm c 

RxV RxV 



We will prove now that ||ofatfc|| s is bounded. To do this, let us apply the Weitzenbock 
formula to the sequence 



Then 



2 

|2\ 



+ \ \V dukWl + Ricc(du k , du k ) =< du k ,dA g u k > g (18) 



3 ^ 2 k ^ 9 ^ + \ \Vdu k \\ 2 g + Ricc(du k ,du k ) + c fc _i| \du k \ \ 2 g + ^9(b k - 1 )(\\du k \\ 2 g )+ 

[ < gradu k , V gradUk b > g +du k (J^j < du k , du k ^ > 9 ]| A =„ fe _ 1 = 

Oc 

< df,du k > g -u k [< dc,du k > g +— < du k ,du k -! > ff ]U=« fc _ 1 (19) 
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where Cfc(x)=c(u fe (x),x) and bfc(x)=b(ufc(x),x). In order to prove that the gradient of the 

sequence is bounded, we evaluate the formula(19) at a maximum point P of ||V«fc||^. 
Because A g (\\u k \\ 2 )(P) > and 6(b k )(\\du k \\ 2 MP) = 0, we get: 



oo 9c 

(co-ero-6o-/3||dufc-i||oo)||dwfc||« < IM/1U+ — — [lMclloo + niax[ — l-H^-iHoo] (20) 

Co V OA 

Inequality (20) implies the following estimates for e small enough, for all k > 1, 

(co - b - r e - PWduk-iWoJWdukWn < fo + 7lM«fc-i||oo (21) 

Elementary properties of homographic recurrent sequences show that | |dtifc_i| 1^ < R{e) 
implies | \du k \ \ (XJ < R(e) for any e E]0,e[ (see 17 for the definition). Let us prove the result 
by induction: since Mo — 0, we obtain from equation (21) that maxy n < - _/°_ £r • 

A simple computation proves that 6 ^°_ < -R(e) and the induction property follows. 

In order to prove that the sequence u k converges uniformly, we will prove the uniform 
convergence of the series w k = u k+ ± — u k . Let us prove that each term of the series{u>fc} 
is bounded in absolute value by the terms of a converging series . Recall that b u (x) = 
b(u(x), x), c u {x) = c(u(x),x) for xGV. Then from 

€A g U n + 0(&u„-i)Wn + c Un _ x u n = f 
eA g u n+l + 6(b Un )u n+ i + c Un u n+ i = /, 

we get: 

eAgW n + O^bu^Jwn + c Un _ Y w n = 6»(6 Mn _ 1 - b Un )u n - (c Mn _ 1 - c Un )u n 

Let: 

g« = °ip Un _ 1 - b Un )u n - (c^ - c Un )u n 

w 2 satisfies: 



n 



-A g w 2 n + 6{b Un )w 2 j2 + c Un w 2 n = w n g n - \ \dw n \\ 2 g 

At a maximum point P of w 2 , w 2 (P) = maxy n w 2 , Aw 2 > and 0(b Un )w 2 = 0. Hence: 
c Un (P)\w n {P)\ < g n {P) and 

( db dc \ 

9n(P)< (maxl^l+maxl^lj \ Wn _ l{P )\ \\du n (P)\\ g 

Since c(A, x) > c > 0. We obtain the estimate for w n : 

i I f r / ,3iin n / v / i <9c / max I/I \ 
max wJ < — ^ (max — — Jit(e) + (max — — ) \ max u> n _i 
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Then: 

max wJ < max ™ n _i 

V c v 



Using the definition of R(e) : 



/3R(e) + 7 = Co - 6 - ero - -A- (22) 

R{e) 

and we obtain for e small enough: 

gg(e) + 7 < c - &o - er _ / i 
c ~~ c Co-R(e) 

This proves that the series w n converges geometrically in C° topology. Finally the se- 
quence u k converges in the same topology and by the elliptic estimate, it will converge in 
C 2 ' a topology (see [7]) and R(e) is a Lipschitz bound for the limit: max| \du e \ \ g < R(e).u e 

is a smooth solution of the equation: 

eA g u e + <b Ue ,gradu e > +c Ue u e = f (23) 

Weitzenbock Formulas In order to obtain regularity results,we need an estimate Vdu £ 
uniform in e. We shall establish a generalization of Weitzenbock identity which we proceed 
to prove. Define an extension of the Laplacian A g to covariant tensor fields u as follows: 

A g u kl ,..k n = — V fc Vfcu;fc lj ..fc n 

A g u ku .. kn = g ij [(V,V^) fci) ikn - rj- (V k u) ki _ kn 
For any C 4 — function u: 

AgVdu = VdA g u + n (Vdu) + K^du) (24) 

where TZq\® 2 T*V ><S> 2 T*V is an endomorphism of the bundle of covariant 2-tensors 

into it self, defined as follows: 

Ko(u) id = R\uj 31 + R 1 ^ + 2g kl R k -- m u ml 

and 72-1 :T*V — > ® 2 T*V, the vector bundle homomorphism: 



■ni 
-jn- 



Rij is the Ricci curvature R^- , R* = g m R n j = Rj 
Then the following two facts are easy to check: (i) TZq maps symmetric tensors into them- 
selves, (ii) the image of 72-iis made up of symmetric tensors. Multiplying equation(24) 
scalarly byVc£u,a simple computation shows that : 



-A g \\Vdu\\ 2 g + \ \V 2 du\\ 2 g =< VdA g u,Vdu > g + < 1Z (Vdu) + TZ 1 (du),Wdu > g . (25) 
The relation (25)is a generalization of the classical relation(18). 
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3.2 Regularity estimates 

Proposition 2 Using the notations of Proposition(l), under assumptions 1-2, there exists 
ae > such that the function ofe, ||Vdw e ||oo, is bounded on any closed subinterval of [0,e[. 

Proof: We take the exterior derivate of equation(23), multiply the result scalarly by 
Vdu £ .We get: and to the solution u e of equation(23) and use this equation : 

e < VdA g u £ ,Vdu e > g + < Vd[(9(b Ue )u e },Vdu £ > g +(V 'd(c Ue u e ) , V 'du £ > g =< Vdf,Vdu £ > g 

To compute < V 'dA g u £ ,\/ ' du £ > 9 ,we apply the relation (25): 

-AJ|Vdw|| 2 ,+£||V 2 dw|| 2 ,+ < Vd[(9(b Ue )u e ], Vdu e > g +(Vd(c Us u e ),Vdu £ > g =< Vdf,Vdu £ > g + 

(26) 

+ < TZ {Vdu) + TZ^du), Vdu > g . 
Let us estimate the third term on the left hand side of equation (26). We have: 

Vd[(9(b Ue )u e ] = 6(b Ue )Vdu £ + B 2 (b Us ;du £ ) + K(b Us ,du £ ) (27) 

where, for any vector field F on V,uj GT*V — >£> 2 (F;cc>) e (g) 2 T*V is the vector bundle 
homomorphism defined by: 

B 2 (F;u) = u(V 2 F(q)), Hue T q V 
and 71: TVxyT*V — ><g> 2 T*V is the homomorphism: 

n (x,u)ij = -K ikj l x k UL 

Note that the images of both homomorphisms are contained in the space of symmetric 
2-tensors. The proof of relation(27) is given in Appendix 2. Multiplying both sides of 
equation(27) scalarly by Vdu £ : 

< Vd[(8(b Ue )u e ],Vdu £ > g =< 6{b Ue )Vdu £ ,Vdu £ > g + < B 2 (b Ue ; du £ )+TZ(b Ue , du £ ),S7du £ > g 

6(b UE )Vdu £ = V bu Vdu £ + i5i(6« £ ; Vdu £ ) 

where for any vector field F on V, B\ is the homomorphism: uo G (S> 2 T*V — >B\(F;uj) G 
<8) 2 T*V defined as follows: 

B 1 {F; u)ij = uj kl V 3 F k + u kj ViF k . 
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We have: 

< 9(b Ue )Vdu £ ,Vdu £ > g =< V hu Vdu £ ,Vdu £ > g + < Bi(b Ue ;Vdu £ ),Vdu £ > g 
< 9(b Us )Vdu £ ,Vdu £ > g = ^9(b u _ :) \\Vdu £ \\ 2 g + < #i(fc Ue ;V 'du £ ) , V 'du £ > g 

< Vd[{9{b Uf )u e },Vdu £ > g = ^e(b Ue) \\Wdu £ \\ 2 g + 

< Bi(b Ue ]Vdu £ ) + B 2 (b U£ ; du £ ) + TZ(b Us ,du £ ), Vdu £ > g 
Finally we get the identity: 

e -k g \\Vdu\\ 2 g + e\\V 2 du\\ 2 g + -9{b Ue) \\Vdu £ \\ 2 g + c Ue \\Vdu £ \\ 2 g + S 1 = S 2 + S 3 

where: 

51 =< B 2 (b Ue ; du £ ) + < Bi(Vdu £ ), Vdu £ > g 

5 2 =< —TZ(b Ue , du £ ) + e[R,o(Vdu t ) + 7Zi(du £ )], Vdu £ > g 

5 3 =< u £ Vdc Ue + dc Ue <S> du £ + du £ ® dc Us , Vdu £ > g + < Vdf, Vdu £ > g 

To prove that 1 1 \/du £ \ \ O0 is bounded we follow the same method as the one used to prove 
the boundedness of Hdu^!^. Let PgV be a point where the functiomxG V->\ \ Vdu £ \ \ 2 (x) 
attains its maximum. Then |A g ||Vdw||^(P) > and \6{b Ue )\\V du £ \\ 2 (P) = 0. 

Hence: 

3 

c Ue (P)\\Vdu £ \\ 2 g (P) <J2\S n (P)\ (28) 

n=l 

To estimate Si, it is convenient to choose an orthonormal coframe field u) l 1 ... 1 u) m in a 
neighborhood of P. Denote by ei,...,e m the corresponding frame field. Then: 

Vdu £ = U af3 u a <g) uj! 3 

where the matrix of functions U a/ 3 is symmetric in a, (3. Moreover, because U a/ g is sym- 
metric we can choose the coframe field uj 1 , ...,uj m so that at P: 

U aP {P) = 0iia^l3 

Vb = B*e a ® c/ 
B l {Vdu £ ) aP = U^Bl + U M Bl 

m 

<B l {Vdu £ ),Vdu £>g = U aP {U ai B} + U Pl Bl) 

a,/3,7=l 
m 

< B 1 (Wdu £ ),Wdu £ > g (P) = U^PW^P^BJ + B^P) 
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< 



B 1 (Vdu £ ),Vdu £ > g (P)| < b J2Uaa(P) 2 = b \\Vdu £ \\ 2 g (P). 



a=l 



To estimate the term B 2 (b Us ,du £ ), for 00 e T*V,X,YeT 9 V,qe V: 



S 2 (6 U£ ,d Me )[X,F] = du £ ((V 2 b) Ue [X,Y]) + du £ Vx 



96 



.cfai £ (Y) + 



du F V 



96 
9A 



.<MX) + dw e ( J^(u e (g),g) ) AVi/.f.Y. V] + 



<9 2 6 



Hence: 

||i3 2 (& Ue ,dw £ 

and 



3 = ||^ £ || 00 .||V 2 6|| 00 + 2||rf M£ ||^ . 

/^II^IU.IIv^iup) 



„db 
V 9A 



<9 2 6 



<9A 2 



+ 



„db 



\Si(P)\ <b \\Vdu E \\ 2 g (P) + R(e)\\Vdu(P)\\ g {\\V 2 b\\ oo + 2R(e 
+ (3\\Vdu £ \\ g (P)} 

I |Si I loo < (b + /3i2(e))||Vdu c ||^ + BR{e) [1 + R{e)] 2 \\Vdu t 



+ R(sf 



d 2 b 



d\ 2 



(29) 



where: 



„db 




<9 2 6 








v Tx 


j 

00 


d\ 2 



} 



P = max{||V 2 6|| 00 , 

Estimate of S 2 (P): 

\S 2 (P)\ < 1^(6)1^ + e(K 2 R(e) + K 3 \\Vdu £ (P)\\ g )]\\Vdu £ (P)\\ g , 

where the constants Ki and K 3 depend only on the curvature of g and K 2 on the covariant 
derivatives of the Ricci curvature tensor. 

Halloo < [K 1 R(e)\\b\\ 00 + e(K 2 R(e) + K 3 \ \ Vdu £ \ U] \ \ Vdu £ \ U (30) 

Estimate of S 3 (P): 

Br 

\S 3 (P)\ < ||Vd«||oo{||V^||oo + C(l + R(e)) 2 + |«(P)^(«(P),P)|.||Vd«|U} 
where C is a constant depending on f and the derivatives of c up to the second order. 

HSjIoo ^ l|V^|| 00 {||Vd/|| 00 + C(l + P(e)) 2 + 7 ||V^|| 00 } (31) 
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The inequalities (28 29 30 31) imply the estimate: 

c HVdUelloo < (6o + PR(e) + 7 + eK 3 ) HV^IU + [BR(s) + C](l + R(e)) 2 + 

Xii?(£)||6||oo + £^ 2 i?(£) + HVd/IU (32) 

The relation(22) implies that: 
/o 



|?f ; . £ [r + ir 3 ]J ||VdM e ||oo < [ J Bi?(e)+C , ](l + i?( £ )) 2 + [K 1 ||6|| 00 + eK 2 ] J R(£) + 

Hence HVdi/eHoo is bounded on any proper subinterval of [0,e[ where e is either e or the 
root of : eR{e)— ^ K if it exists. 



4 Study of limit when the viscosity parameter con- 
verges to zero 

We will prove here that the first order partial differential equation (16) has a Lipschitz 
continuous solution obtained as a limit of the sequence u e when e converges to zero. 

Theorem 2 1. Under the assumptions 1-2: (i) Any sequence {u £n \ e n ->0 as nG N, 
tends to 00} contains a sub-sequence that converges to a solution of equation (16) 
of class C x in the C 1 — topology. Hence there exists C 1 solutions to the first order 
partial differential equation (16). The solutions of the equation (16)which are limits 
of sequences {u £n \ e n ->0 as nG N, tends to 00} are called viscosity solutions. 

2. Any viscosity solution has Lipschitz continuous derivatives with a Lipschitz bound 
independent of the solution. 

3. For any M >0, there exists a neighborhood N (M) of in C 0,1 (V) such that for any 
/G Af(M),the equation(16) has at most one Lipschitz continuous solution, having M 
as Lipschitz bound. 

Proof : (i)To study the limit of u e when e goes to zero, we use the same type of 

estimates as previously. For all e > 0, sup| \du £ \ \ g < R(e). It is easy to see that R:[0,e] 

v 

— >M + is a decreasing function. So any sequence {u e J e n j as n | 00} contains a 

subsequence converging in the C°-topology to a Lipschitz continuous function u such 

that sup||du|| g < R(0). Any limit u of u e , as e tends to 0, satisfies the weak equation 
v 

where the integrals are taken with respect to the volume measure associated to the metric 
g: 

V0eifi(F), / [u{6{b u )<t>) + c u u<t>\ = [ f<j> 
Jv Jv 
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and almost everywhere on V, the equation: 

du(x) .b(u(x) , x) + c(u(x), x)u(x) = f(x), x G V 

Because Vdu £ remains bounded as e goes to zero, we can find a subsequence of {du £n }, 
that converges in the sup norm. Then a second subsequence {u £n } will converge in the 
CHopology to a (^solution of equation (16). 
(ii) Follows from the preceding. If u is a viscosity solution, then: 

Lip(u) < 

Jo 



As in the linear case we have the formula: 

,o 

u(x) — / K(t,x)dt, x G V, where (33) 

J — oo 

K(t,x) = f(v?"(t, a;)) exp[- c(w(y9 M (s, x)), V u (s, x))ds] 

Here {</?"| tG R} is the flow of the lipschitz continuous field b u . 

Now we shall prove statement (ii). Assume that f belongs to the neighborhood JF(M) 
and that equation( 16) has two distinct solutions u ,Ui both lipschitz continuous with 
lipschitz constant M. Let v=ui-u . Then using the relation (33) we get: 

v(x) = f {[f( V ul (t,x) - f(cp^(t,x)]e-^< u ^ ul ^ x ^ ul ^ ds + 

J — oo 

f((p uo (t : x))[e~ f° c ( ul (' pUl ( s ' x ^'' pUl ( s ' x ^ ds — e — c ^ ti0 ( S ' X ))'V"° (s,aj)))ds j ^34^ 

where <p Ui denotes the flow of the Lipschitz continuous field b Ui , i=0,l. To proceed 
we need a careful estimate of the max d g (ip U2 (t, x), <p ui (t, x)).d g is the distance function 

on V defined by the Riemannian metric g. Let u r =(l-r)ui+ru2 , rG [0,1] and denote 
by {ipl\t G R} the flow of b Ur (instead of the more complicated notation <^" r ). Assume 
temporarily that uo,ui and hence u r are smooth. Denote by ?/>:Rx Vx [0,1] — >TV the 
mapping:^ (t,x,r) = -^y^. tp satisfies the variations equation: 

V t tp(t,x,r) = {V^ t ,x,r)b){u r ((p r (t ) x)),(p r {t,x)) + [v(<p r (t,x) + du r (ip r (t, x)).ip(t, x, r)] 

— (u r (<p r (t,x)),<p r (t,x))) (35) 



Multiplying scalarly both sides of equation (35) by ip(t,x,r), we get: 

< V t ip(t,x,r),ip(t,x,r) > g = (V v>( (i i:Cir )6)(M r (v? r (t,x)), <p r (t,x)),i(>(t,x,r) > g + 
db 

< — (u r (ip r (t, x)), <p r (t, x)),if)(t, x, r) > g .[v(ip r (t, x)) + du r ((p r (t, x)).ip(t, x, r)\ (36) 
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Equation(36) in turn implies the inequality: 

~ll^ie<6d||^li; + )9[IHIoo + M||^|U||^|| fl . (37) 

because u admits M as Lipschitz constant. Gronwall's lemma applied to the inequality(37) 
gives : 

e t(bo+f3M) 

I m, x,r)\\ g < P\\v\ \oo7——TTj for all (t, x, r)mR + xVx [0, 1] 

and ■?/'(0,x,r) = ^ !i =0, because (p r is the identity map of V. For negative times one uses a 
time reversal and changes b into -b. Finally we get: 

e \t\(b +/3M) 

I Mi, x, r)\\ g < 0\ Moo bo + (3M for all (f , x, r) in R x V x [0, 1]. 
Because re[0,l] — >ip(t,x,r) e V, is a path joining <p u °(t, x) to (p Ul (t,x): 

/•l |t|(6 +M/3) 

d ff (^(t,x),^(t,x))< J \m,x,r)\\ g dr<(3\\v\\ 00 bQ + f3M (38) 

In the case where u ,Ui are just Lipschitz continuous with Lipschitz bound M, we apply 
lemma(l) to approximate Uj,i=0,l, by a sequence of smooth functions {uj 5 „| nG N}, Uj jn 
with Lipschitz bound and take the limit in the inequality (38) corresponding to 

Ui >n ,i=0,l. 

The equation (34) gives us the following estimate: 

rO e \t\(b+(3M)-co\t\ ( dc \ 

\\v\U<\\v\U{(5Li P {f) J ^ bQ + (3M ^ + /3||/|U^Mmax|-|+max||&|| g j. 

fO rO \s\(b Q +M(3) a rO 

IsL wTw dsdt+ L tetc " dt} 

This implies that: 

(3 ( dc \ 1 dc 

I If I loo < \\v Hoo-llf lie* - ih ; t-tt; ( 1 + Mmax|— - 1 + max||rfc|L H — ^max|-— 11 

We chose as neighborhood J\f(M) the set : 

3 ( dc \ 1 dc 

{g G C^{V)\ \\g\\ c ^ < 1/t , 1 + Mmax| — | + max||dc|L +^max| — |} 

We conclude with some remarks. 

Remark 1: The equation(16) can be interpreted as follows. Let X denote the vector 
field onlxV whose first and second components at (A,x)g IxV are respectively f(x)- 
c(A,x)A and b(A,x). Then a function u:V — >R is a solution of equation(16)if and only if 
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the graph of u is an invariant manifold for X . Hence one can apply the results of [8, 4, 15] 
to show the existence of a solution u. The advantage of our method is that it gives precise 
estimates insuring the existence of a solution. 

Remark 2: Our result about the existence of a solution for equation( 16) states that 
the perturbation" < b u ,gradu > +c u — f " of the equation "< b u ,gradu > +c u = 0" 
has a solution near the obvious solution of "< b u ,gradu > +c u = 0" provided that the 
perturbation f is small enough. Hence our results (as well as those of the papers quoted 
in Remark 1) are purely local. 

Remark 3: An interesting question albeit a difficult one is to study the structure of 
the set of solutions of equation(16), in particular to find the number of solutions if the set 
is discrete. Indeed, in the linear case, there is a unique solution and the previous theorem 
2 proves that there exists a unique regular solution in a small neighborhood of the zero 
solution. Outside, it is not clear how multiple solutions will appear. Is it possible to give 
conditions on the recurrent sets of the vector field that would imply a finite or countable 
number of solutions u. 

Remark 4: It is easy to adapt our methods to the obtain solutions of the equations: 

eA g u t (x)+ < b(u £ (x) , x) , gradu £ (x) > g +c(u € (x),x) = 0,x e V 

< b(u(x),x), gradu(x) > g +c(u(x),x) = 0, x G V (39) 



provided that an approximate solution w:V — >R, of equation(39) is known, such that the 

inf 

v 



error is small enough and inf dc W*M i s sufficiently big. 



5 Some examples 

We will give some explicit examples that are usefull to understand how various can be 
the solutions of first order partial differential equations. 

Example 1: Consider real numbers K,a, (3, K large and, (3, a small and a positive function 
/ on the 1-dimensional torus S 1 . Define the function c and the vector field b as follows: 

c(u, x) = K + K/3(l + a cos(x))u/2 — j3au/2 shire 

d 

b(u, x) — (1 + P(l + au cos(:r))— — 

ox 



The equation : 
is equivalent to 
where 



Ou 

b(u(x),x) — (x) + c(u(x),x)u(x) = f(x) , x e V (40) 
ox 

dU 

^KV = J (41) 

u 2 

U = u + (3{l + acos(x))— (42) 
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Because f > 0, there exists a positive periodic solution of equation ( 41). If K is big , 
so is min c. Then equation (42) has exactly two real regular solutions, if the quantity: 

[0,2tt] 

A(x) = 1 + 2/5(1 + acos(x))U is positive. It is the case if (3 is small enough. In this 
case, one solution is positive, the other is negative. But on the other hand, when A(x) 
vanishes, the solutions are not differentiable any more. When A(x) changes sign, the 
solutions becomes complex and in the real domain, there are pieces of intervals with no 
solutions. 

In the context of this paper, how this example fit with proposition-1? Since the 
solution u is bounded a priori , the assumptions 1-2 are satisfied. In this context, one of 
the two solutions is selected by the recurent process. 

For the singular perturbation problem associated to equation(23), what is the solution 
of the first order partial differential equation, that will be selected at the limit, when the 
small parameter tends to zero: 

By the maximum principle if the first term uq of the recurrent sequence Uk is positive, 
all the functions in the sequence are positive and so is any limit of u e as e goes to zero. 
Hence the sequence converge to a positive solution. Can one find an initial data, such that 
the sequence converge to the other solution. How to describe the basin of attraction 
of such solution associated to the map k — >k + 1. 

Example 2: In this second example, conditions 1-2 of theorem (2) are satisfies and 
the quantities c and b are now bound with the variable u. Consider indeed, 

e~ u2 - 1 e -- " 2 - 1 

c(x, u) — K(l — (3(1 + a cos(x))) (5a 12 sin a; 

u u 

2 d 

b(x,u) = (1 + (3(1 + acos(x))ue~ u )—. 

ox 

Equation 40 corresponding is equivalent, as in the previous example to the equation 41 
with 

e~" 2 - 1 

[/ = u- (3/2(1 + a cos x) . (43) 

u 

It is easy to check that for K large, c is big. The conditions 1-2 are satisfied, if (3 

is chosen small enough, because the function g : u — > and its first derivative are 

bounded. In that case there exists generically several solutions to equation (43). 

Remark: This last example shows how the coefficients of the partial differential 
equation can control the behavior of the sequences u e as e tends to zero. There exists 
other types of examples which exhibit shock phenomena: this is given for example in 
dimension 1, with Burgers equation, with a forcing term (see [12]). 

As mentioned above, all solutions of first order partial differential equations are not 
equivalent with respect to the singular perturbation process and in general there are no 
criteria to distinguish among them. To our knowledge, there is no theory about the 
stability of such solutions. We shall not address these questions here. 
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5.1 Ergodic fields 

In this section, we analyze very briefly the case where the field b is ergodic for the 
Riemannian measure, which implies that divb=0 . In this case, we obtain for the solution 
of the first order P.D.E. (9), we prove an averaging result without assuming that u is 
differentiable. Remark that u cannot be constant if the ratio - is not constant. 

c 

The important question is what determine the value of the solution: is it enough to 
give a value of the solution u at an arbitrary point, and due to the ergodicity, the field will 
spread the information ? Due to the recurrence of the field, we expect that it should be 
the case, because along the characteristic, which dense trajectories, the value at a specific 
point will propagate on the manifold. This situation suggest also to study solutions of 
equation (9), when the field b is a Hamiltonian system. The phenomenon involves here 
should be really different from the case of hyperbolic fields. 

No limits are expected when t tends to infinity, since all trajectories are dense in the 

min|/| max|/| 

manifold : the solution u will oscillate between — — and —. — . Using the ergodicity, we 



max c mm c 

V V 



can compute lim^ +00 | J"' u{x{t))dt explicitly. 

For simplicity, we suppose that c is a constant function. Let PgV and let x(t) be the 
trajectory passing through P at time 0. Then: 

u{x{t)) = u (P)e-ti< x ^ ds + f f{x(r))e-^ c ^ ds dT 

Jo 

Averaging u, we get 

- [ u(x(t))dt = -u(x(t )) [ e-< s - u ^du + - [ [ f{x(T))e-< s - T ^dTds 
t Jo t Jo t Jo Jo 

and changing the order of integration, 

- / / f{x(T))e- c ^dTds = - f f(x(r)) [ e-< s ^dsdr 
t Jo Jo t J J t 

= i£Z(£M> (1 _ e -^> )Al 

The last integral behaves like| J * dr and we conclude by using the ergodic properties 
of b: 

1 f* iff 
lim - / u(x(t))dt = UTr , V 
t-^+oo t J y w; Vol(V) c 

6 When the zero order term vanishes 

In this paragraph, we emphasize the importance of the positivity of the lower order term of 
the first order PDE studied in the previous section. In particular, by considering a simple 



25 



equation, we will prove that the results proved in the previous section like theorem-2 are 
no longer true. Consider the linear partial differential operator on V: 



L £ = eA g + 9{b) 

and the associated equation: 

eA g u £ + e{b)u £ = /, (44) 

f is given function on V. Such an equation has a solution u e if and only if f is orthogonal 
to the kernel of the adjoint L* = eA g — 9(b) — div g b of L £ . On the other hand if there is 
a solution, it is not unique since the kernel of L e is not reduced to zero. All the constant 
functions, for example, are in this kernel. In order to avoid these difficult questions we 
shall replace the equation (44)by the following one: 

eA g u £ + 6(b)u £ + eu £ = f (45) 

Equation(45) has a unique solution u £ for e > 0. We shall prove that there exists choices 
of b and / such that the sequence u e does not remain bounded as e tends to 0. As 
expected, the first order partial differential equation (9) has no bounded solutions u for a 
generic f. 

Proposition 3 Suppose that a is a critical point of the field b and that:l)either the lin- 
ear part B of the field b at a is symmetric positive with respect to the Euclidean metric 
g a induced by g on T a V, 2) or more generally B=S+A, S,A<EEnd(T a V), such that S is 
symmetric positive definite with respect to g a and SA is antisymmetric with respect to g a . 
Assume also that f(a) ^ 0. Then u e does not stay bounded when e goes to zero. More 
precisely, it diverges to infinity at the point a. 

Proof:To begin with, we need some remarks in case 2). The conditions on B imply: 
(i)The function: xGT a V — ><Sx,x> 9a is invariant under the action of the one parameter 
group e tA generated by A.(ii) The trace of A, tr s A with respect to g a is zero. In fact if 
A* denotes the adjoint of A with respect to the scalar product < , > ga , SA+A*S=0. So 
SAS~ 1 +A*=0 and tr fl (SAS~ 1 )+tr 9 A*=0. But this gives: 2tr 9 A=0. 

The proof proceeds by contradiction : suppose that the sequence u e is bounded. 
Consider a normal geodesic coordinate system at point :B— >R, x*(a)=0, 

1< % < m, the image of which is the ball B m (0,r), with center and radius r in R m . The 
rescaled function: 

v e (x) = u e (y/ex)) 

defined on B m (0, -^), satisfies the following equation: 

A 9e v e (x) + 9( b -^P-)v e (x) + ev e (x) = f(Jix) , x e B m (0, 

g e denotes the rescaled metric. The coefficients of this equation are bounded. When 
e tends to 0, the rescaled metric g e converges to canonical metric on R m and the field 



26 



converges to its linear part at 0, i.e. the linear field YlTj=i ,j% 3 -£j , Bjj=^-(0), 
uniformly on every compact in R m . For any sequence {e n \ nG N} converging to 0, it is 
possible to extract a subsequence still denoted by {e n \ nG N} for simplicity, so that {v €n } 
converges also uniformly on every compact in R m to a bounded solution v of the following 
partial differential equation(see [7]): 

m ft 

A e v + J2 r> >-r> J -^ = /(°) on R m (46) 

Here A e = - £J 

=i (6>x«)2 ■ us now introduce the potential function: U:T a V >R: 



C/(x) =< > 9a = -^B 
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Let ^ = exp Equation 46 becomes after some computations 

A e (vijj) - vA e ^ = /(0)V> on M m (47) 

This can also be written as: 

A e (vip) = vA e ^j + f(0)ij; on R m 

The right hand-side of this equation is a continuous function on R m tending to zero 
exponentially at oo. Because the function vip is bounded on R m , its partial derivatives 
are bounded on R m (see [7]). 

Finally consider the change of function w=v0- Then equation(47) becomes: 

^A e w - wA^ = /(0)V 2 

To find a contradiction, we will integrate the last equation on a large ball B(R) of center 
and radius R in R m with respect to Lebesgue measure. Using the Green formula on 
equation the previous equation, we obtain calling da R the area measure on the sphere 
dB(R) and n the unit exterior normal: 

f (^ - w^)da R = [ fm 2 on R m 
JdB(R) on on J B(R) 

Since w and dw are bounded, and ip is a Gaussian , the integral j dB ^{^^ ~ w ^)do'R 
converges to zero as R tends to infinity but the integral f Rm ip 2 exists and is not zero. 
Hence we arrived at a contradiction if /(0) 7^ 0. So the sequence u e is not bounded and 
no bounded solution for the first order PDE is obtained by this method. 

Let us prove now the proposition in the second case. Then v is solution of 

Av+ H S '-r' J ~^ + E - V -r> J -^ = /(0) + on IT 

i,j=l i,j=l 
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Consider now the potential U:T a V — >W , associated to the endomorphism S, given by : 

U(x) =< Sx,x> g = J2^- 

and let ^ = exp = ^-. Using the transformation v — >w=vip as before, we see that equation 
transforms into 

m P) 

4>Aw - wA4> = /(0)^ 2 - *p 2 A iJ xi JTl 011 Rm ( 48 ) 

We follow the same steps as previously, the only difference is the existence of the second 
term on the right hand-side of equation 48. Integrate both sides of the last equation on 
the solid ellipsoid E(R)={x| U(x)<R 2 }, with respect to Lebesgue measure: 

I (4>Aw - wA^) = [ f(0)^ 2 - [ ^ 2 J2 Aijxtp, 
Je{r) Je(r) Je(r) f-j± 1 ox 

Using Green's formula to estimate the term on the left hand-side, we get: 

/ ^Aw-wAij)ds dE(R) = [ fm 2 ~ [ ^it^TTi 
Je(r) Je(r) Je(r) -~Z x ax 

where ds e(r) denotes the area measure on the ellipsoid dE(R) = {x\ U(x)=R 2 }. The 
second integral on the right hand-side of equation is zero. Assuming this for the moment 
we get a contradiction in exactly the same way as in the case when B is symmetric. 

To see that the second integral on the right hand-side is zero, denote by A the linear 
vector field on W 71 having Y^j=i ^i,j x ^ as the i th component. The condition on the matrix 

A implies that: 9{A)U = 0. This implies that E(R) is invariant for the flow of A. Hence: 

0= f eiAj^vdx 1 A .... A dx m ] = [ { [(0( ~A)v)i) 2 + v OiAjl^dx 1 A .... A dx m 

+ ^ 2 v eiA^dx 1 A .... Adx m ]} 

But: 6(Aj(ilj 2 ) = -i/j 2 6(A)U = Oand G( A) [dx 1 A . . . . Adx m ] = trA^x 1 A....Adx m ] = O.And 
0(A)v - , A u x^. 

Remarks: In both cases 1) and 2) a is a hyperbolic repealer for the flow of b. In fact 
the function U is a Liapunov function for the flow of b at a in both cases. 

Generically, when f does not vanishes with b, equation 9 does not have bound solutions. 
Using the characteristics solution and equation 14, any solution if it exists, should blow 
up near 0, a zero of b like "Clog||x||, where \ \x\ \ is the distance from the point x to 0. 
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6.1 The gradient case with symmetries 



Consider a gradient vector field b = grad(f), then equation 

eA g u e + 9(b)u e = f 

can be transformed into 

div((f> 2 gradu e ) = f(p 2 
& = 

The integral condition imposes on a compact manifold that J v f<p 2 = 0. Suppose that the 
manifold is a sphere, f is odd and <fi is even. The integral condition is then satisfied by 
the symmetries assumptions. Assume that <fi has a minimum where / is not zero. We can 
apply locally the results of Proposition-3 to show that the sequence u e is not bounded. 

6.2 Negative case 

We finish this section by mentioning an open problem concerning the construction of a 
sequence of approximate solution of the following equation: 

A g u t + < b, Vu e > —e~ 2 u t = f on V 

This problem is related to the linear wave equation on manifold when b = 0. It is known 
that there exists a subsequence of u e that concentrates along a stable close geodesic (see 
Babich-Lazutkin, [1]) of the manifold. We ask the following question: when e tends to 
zero, is it possible to construct a subsequence that concentrate along the limit sets of the 
vector field, instead of the geodesies? moreover what are the minimal assumptions on the 
vector field b that insure the concentration of the sequence of eigenfunctions ? 

7 Conclusion 

We have studied the scalar case. But our methods could be extended to study systems of 
the form: 

TO 

eAgU + y^ j B n (u(x),x)— + C(u(x),x)u(x) = f(x),xE V 

n=l 

u ( x } ,x }~fa™ C(u(x),x)u(x) — f(x),x G V 

n=l 

where u,f:V — >R JV and B n ,C :RxV — >End(lR Ar ) are given matrix functions. C must 
satisfy positivity conditions (see [4], [5]). 
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8 Appendix 1 



In this appendix we give a proof of lemma(l). We embed (V,g) isometrically into a 
Euclidean space R n endowed with its canonical metric. There exists a r* >0, depending 
only on the normal curvatures of V in M n , such that for any rG]0,r*[, V has a tubular 
neighborhood (T r ,7r r ,V) of V in MJ 1 of the following type:for any xGV, 7r r (x) is the open 
ball of center x and radius r in the affine space N X V, passing through x and orthogonal to 
V at x. It is easy to see that there exist a function c:[0,r*[ — >[0,+oo[, such that c(r) tends 
to as r tends to and max{T7r r (v) | v tangent to T r and ||v||r"=1}< l+c(r). 

Call k r : T r — >R, the function ko7r r . As a function on the manifold T r , k r is lipshitz 
continuous with lipschitz bound M(l+c(r)). This means that for any couple x,y in T r , for 
anyC 1 curve7 : [0, 1]— >T r , such that t(0)=x, 7(l)=y, | |fc r (ic) -k r {y)\ k™ < M(l+c(r))x 
length of 7. Rademacher's theorem implies that k r belongs to the Sobolev space W 1,0O (T r ). 

Choose an r such that 2r<r* and that Mc(2r)<e. Let a:M n — >[0,+oo[ be a C°° function 
with support in the ball of center and radius r in R n and such that J Rn a(x)dx=l. 
Consider the convolution h r :T r — >W: 

h r (x) = / k r (y)a (x — y) dy = / k r (x — z)a(z)dz 

Clearly h r is C°°. If x,yGT r and ||x-y||]Rn<r, the straight line segment joining x-z to y-z is 
contained in T 2r , for any z such that | |z| |r« <r. Then for all these z — z) — k r (y — 

z) ||r« < (M + e)\\x - y\\mn. Hence \\h r (x)— h r (y)\\mp < (M + e)\\x — y\\mn. This implies 
that \\dh r \ \ <M + e. 

If xG T r and z is such that ||z||]Rn<r, the segment joining x and x-z is contained in 
T 2r and \\k r (x — z) — k r {x)\\^ < (M + e)r. Hence: 

\\h r (x) - k r {x)\\Rn < (M + e)r for all x G T r , 

By taking r sufficiently small the restriction h of h r to V will belong to il and will have 
M+e as Lipschitz bound. 

9 Appendix 2 

In this second appendix, we give the computations related to formula 27. Let u be a 
1-form. We have: 

6(b)u = V b cu + A b uj 

where: 

A b iu[X] = cu[Vxb] 

for any vector field X. Hence: 

V9{b)u = VVju + VA b u 
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We compute both terms on the right hand side using a coordinate system (x 1 ,...^™). 

VVbCj = WiWbUjdx 1 <g> dx j 
ViV b u = V b ViU + ViVVjU + R{ei, b)uj 
R(e h b)u = -R'^ l b k uJidx j 
ViV h Uj = VbViLUj + V t b k V k oo 3 - Rr k ' l b k cui 
VV 6 cj = VbVto + VVVjW - RgjWwdx* <8> dx j 

A b uo = uJiVjb l dx j . So it is the contraction of the 1-covariant tensor uo and the (1,1) 
variant tensor Vfrand: 

V l (A b u) j = ViU k Vjb k + uj k ViVjb k 

Finally: 

V i e{b)u j = V b V z uj j + V t u k Vjb k + Vib k V k Uj - R' ik ] l b k uji + uj k V t Vjb k 
But we have, if a — a^dx 1 <S> dx^ is a 2- covariant tensor: 

db k db k 
(9(b)a) lj = e(b)oij + a kj — + a lk — 

(0(6)a)ij = (V 6 a)ij + a kj Vib k + a ik Vjb k 
Applying this formula to a = Vcj : 

(0(&)Vw)ij = VbViWj + ViU k Vjb k + V k UjVib k 

We get: 

V0(&)w = #(&) Vcj + w) + 5(6; w) 

where: 

nib-u^^-RHfb^ 

B{b;u) ij =u k V i Vjb k = u{V%) 
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